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equation of state for non-equilibrium systems belonging to the universality class of 
directed percolation (Gribov process) to second order in e = 4 — d. By introducing a 
parametric representation the result can be written to this order in a very simple form. 
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1. Introduction 

Numerous physical systems that show continuous non-equilibrium phase transitions 
display near the transition point space-time structures that can be described by 
directed percolation clusters. Well-known examples are epidemic processes without 
immunization poisoning of catalytic surfaces H, roughening transitions in growth 
processes |§, Schlogl's reaction models j|, f|, contact processes || |7j and certain cellular 
automata || |§. The model known in particle physics as reggeon field theory JK| [□], [L^ 



has also been related to directed percolation and epidemic processes |13], [TJ], [15], [L6| . All 
these systems share the property that their critical behaviour is characterized by the 
existence of an absorbing state. This property determines their universality class || [17] . 

The behaviour of these systems near the transition point is characterized by universal 
quantities such as critical exponents. The exponents for directed percolation have been 
calculated to second order in e = 4 — d and — numerically — for d — 1 and d = 2 (for 
recent results see, e.g., fi"8| , [I~9| , p0[). It is well-known from the theory of equilibrium 
critical phenomena that universality not only holds for critical exponents but also for 
certain amplitude ratios and scaling functions |21], 22, |2J|. An example for a universal 



relationship between physical quantities near a critical point is the equation of state. In 
the present paper we derive an equation that describes the stationary state of epidemic 
processes with an absorbing state near the transition point. 

Before we turn to the mesoscopic description used in the analysis presented below 
let us briefly discuss two examples for microscopic models that belong to the universality 
class of directed percolation. Schlogl's autocatalytic reaction scheme @] is defined by 
the reaction equations 



X + A 2X, X B (1) 

k[ k 3 

where the concentration of A's and B's is kept constant (e.g., by reservoirs). Here two 
remarks are in order: (i) If A — B (i.e., the particles are of the same type) and their 
concentration is not kept constant (but governed by the reaction-diffusion dynamics) the 
system belongs to a different universality class. In that case the total particle density 



is conserved and affects the long-time dynamics of the system [^J, pq] . (ii) Only if the 
reaction B — > X has zero probability the system has an absorbing state. In the case 
of chemical systems it is natural to take this process into account since back reactions 
can in general not be avoided. Thus the B-particles represent a source for X. Below a 
critical value of the reaction rate &2 the density of X-particles is nonvanishing in the 
stationary state for k 3 — > 0. 

In the theory of infectious diseases |]J one considers two types of individuals: 
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Susceptibles S who can catch the disease and infectives I who have the disease and 
can transmit it. The dynamics of the system is characterized by the following processes: 

(i) Healthy (susceptible) individuals may catch the disease by contact with infectives. 
Infected individuals may recover spontaneously and thereby become susceptible 
again. 

(ii) Interaction of infected individuals can lead to recovery ("saturation"), 
(hi) Diffusive spreading of the disease. 

(iv) Susceptibles may catch the disease spontaneously. 

The process (iv) models for instance a continuos flow of germs into the system. Note that 
the system has an absorbing state if (and only if) the process (iv) has zero probability. 

A mesoscopic description for the dynamics near the transition point requires a single 
density n(r, t) which represents the slow degrees of freedom. In the case of Schlogl's 
model ([I]) n corresponds to the density of X-particles, in the context of poisoning of a 
catalytic surface n is the density of vacant sites and in epidemic processes it represents 
the density of infectives. The density satisfies the Langevin equation 

d t n = AV 2 n + R[n\n + \q + ( (2) 

where the reaction rate R[n] can be expanded in powers of n. In order to investigate 
the universal properties near the transition point we only need the two leading terms, 

R[n] = A(r + §n) (3) 

where g > and r is a temperature-like critical parameter. For q = the system has 
an absorbing state with n — if the random force ( is multiplicative: 

(C(r, t)C(r', 0) = A^n(r, t)5(v - r')5(t - t') (4) 

where we have again neglected irrelevant higher powers in n. In the case of the reaction- 
diffusion system (|]) the (constant) source term q vanishes for a zero reaction rate k%. 
By a simple rescaling it is possible to render the coupling constants g and g equal. 



2. Renormalized field theory and perturbational calculation of the 
equation of state 



Hereafter we work with the dynamic functional 0, |2P|, |5D|, |3T], ^2 



J[s, s] 



dt I d d rs 



d t s + A(r - V 2 )s + — s(s - s) - Xq 



(5) 



where s oc n and s denotes a Martin-Siggia-Rose response field [Q. Response and 
correlation functions can be computed by integrating products of the fields s and 
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s with the weight exp(— J). We are especially interested in the average M = (s) 
which corresponds to the mean particle density. The external field Xq describes the 
spontaneous creation of particles and breaks the symmetry s(t) <-> —s(—t). 

In order to obtain the equation of state — i.e., M as a function of r and q in the 
stationary state — we perform the shift s = M + <fr und determine M by the no-tadpole 
requirement (0) = 0. Since the stationary state is homogeneous the shift leads to the 
functional 

JM <f>\ = Jo$, 0] + Jg[4>, 0] + J$, 0] (6) 

where 6 — s, 



J Q [(f>, 0] = / dt / d d r\<f)[M{TM + |M) - q\ (7) 



and 



Jg[$, 



Ji 



dt I d d r^)[d t (t) + A(r + gM — V 2 )0 - yM0 



dt / dV 



The equation of state can now be written in the form 
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Xq = XM(t + -M) - 




(8) 



(9) 



(10) 



where the bubble represents the sum of all one-particle irreducible Feynman diagrams 
of the field theory with one external 0-leg. Details of the diagrammatic analysis 

are given in the appendix. After dimensional regularization of ultraviolet divergences 
one obtains at two-loop order 

G e g 2 T l ~ e/2 



q 



M 



g l r-- 



gM—^h + (2h 



e(2 - e) 24 
where G e = T(l + e/2)/(47r) ci / 2 is a geometrical factor and 
f = r + gM. 



3/ 2 r 



(11) 



(12) 



The integrals I\ and I2 are defined in the appendix. 

The poles in e can be absorbed into renormalizations of the fields and the coupling 
constants. The renormalizations read 



0^0 = Z 1/2 (p (j) -> . 
q - q = Z?ZWq X 
g^g = Z- l ' 2 Z^Z g g 



= Z 1/2 (p M 

* x = z- l z x x 

G e g 2 = ufi e 



► M = Z 1/2 M 

T — > f = Z^ l Z T T 



(13) 
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where /i denotes an external momentum scale. The required renormalization factors 
have been determined to two loop order in Ref. ||. The Z-factors are 

ry U f 7 n 9 l 4 \ ^ 

Z = H h 3 + - In - 

4e Ve 2 3/ 326 

, u /13 31 35, 4\ w 2 

Z A = 1 + — + + — In - (14) 

8e V e 4 2 3/ 128e V ; 

u /I 5 \ , „ it / 5 7 \ 2 

^r = l + 7^+hr^-7T7^U Z q = l + -+[ — - — )u 2 . 



2e V2e 2 32e/ 9 e V4e 2 16e/ 



3. Scaling form of the equation of state 



In order to improve the perturbative result (JTTJ) by the renormalization group we exploit 
as usual the invariance of the bare field theory with respect the a variation of \i at fixed 
bare parameters, i.e., 



q(r,M,g) = 0. 



(15) 



This equation — expressed by renormalized quantities — is the renormalization group 
equation 



Hdp + kt8 t - ]pMd M + (3d u + ( - -7] q(r, M, u; //) 
with the Wilson functions (to two-loop order) 







(16) 



k(u) 

C(«) 

/?(«) 



dlnZ 



dln/i 
d In t 



u 
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u 



dln/i 
dlnA 



-- + 3 --In- — 



3m 



3/ 16 
49 35, 4\ u 2 
T + Y ln 3j64 



dln/i 

du 



u 



+ 



d In /i 



u 



17 , 4 
In - 

2 3 



3u 



it 

128 



(17) 



4\ u 

e + — - (169 + 106 In- 



3/ 128 



The renormalization group equation ( [IB] ) can be solved by characteristics. At the fixed 
point (w* with (3{u+) = 0) 



2e 
"3 



/169 53, 4\2e , 2 

+ — In- — + 0(e 

V192 96 3/ 3 v 



;is) 



the result (combined with dimensional analysis) reads 

q(r, M, u; fi) = ^H^^q^r 1 '^ ^ d/2 r {d+ ^ /2 M, u*; l) (19) 
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with the critical exponents 
V = 7(«*) = - 



25 161 , 4\ . 

1 n- e + 0(e 2 

288 144 3/ V 1 



1 _ 1 e 
2 - k(m*) ~ 2 + 16 



2 = 2 + CK) = 2 - 
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1 + 
1 + 



/107 17 4 

V288 ~ 144 11 3 

67 59 , 4 
1 In - 

288 144 3 



(20) 





5 = 



i-i 

6 



1 



11 53 , 4\ . 2 , 

288 - 144 ln 3) e + ° (e) 



1+^-^=2 + 1 

/3 3 



85 



+ 



53 



ln- 

288 144 3 



e + 0(e 2 ) 



and the scaling function 



g(r,M,^,/i) = M 



6 



r + - -f ( A - 5 ln(f /^ 2 ) - — (ln(f //z 2 )) 2 



12 



+ ^<?M (l + / - 2 ln(f /f) + (ln(f //i 2 )) ; 



where 



A = 1 + e 



13 53 , 4 1 / 
+ In In 3 - - 

288 144 3 4 3 



B = 1 + e 1 In- . 

,288 144 3/ 



(21) 

(22) 
(23) 



An integral representation for / = —2.3439072 ... is given in the appendix. In order to 
simplify the writing we will use hereafter the dimensionless quantities defined by 



/i~ 2 r -> r H~ 2 gM -> M 



H gq^q. 



(24) 



While equation (21) cannot be used directly to investigate the equation of state in 
the critical region (r, M — » 0) one may map the critical region by an appropriate choice 
of the flow parameter I in (|i~9"D on scales on which perturbation theory can be applied. 
For I = M V IP equations (21) and (|I~9"D yield the scaling form 

aq = M s F(brM- 1/p ) (25) 

with the scaling variables x = brM~ x ^ and z = aqM~ s and the universal scaling 
function 

F(x) = x + 1 + %K \(x + 2) ln(x + 2) - 2(x + 1) In 2] 
6 



+ 72 



[x + 4)(ln(x + 2)) 2 - 4(x + l)(ln2) 2 + 0(e 3 ) 



where 



r , , / 85 29 
K = 1 + e I 



- , « 53 

, H In 2 

V288 72 144 



In 3 



1 + e 0.1699728 



(26) 
(27) 



The coefficients a and b in equation (|25]) are 

^2 



a = 2 
6 = 2 



1 + 
1 + 



e e 
3 

1 + ln2 
6 



85 53 89 
18 V48 + 12 11 ~ 24 



ln3--/) + 0(e 3 



(2? 



13 89 



e + — In 3 - 3/ 

36 V48 24 



,99 59 _ _ 53 , * _ _ 

+ ( 1 In 2 In 3) In 2 

v 16 12 24 ' 



O (e a 



(29) 
(30) 



They have been introduced in order to normalize the scaling function, i.e., 

F(0) = 1 F(-1) = 0. 

The equation of state (p5|) satisfies the correct scaling behaviour in the critical 
region but the scaling function (26) can only be used if x is not too large. Since q is an 
analytic function of M for positive r (for equilibrium systems this was first discussed 
by Griffiths [Qj) the approximation (26) breaks down for M <C t@. In order to study 
this limit one may choose the flow parameter in (|T9| ) as I = r u . One arrives at 

aq = (brY & F(M(bT)-P) (31) 

with the scaling variables y = M(6r) _/3 = and z = aqi^br)^ 136 = zy s and the scaling 
function 



F(y) = y \ 1 + y + -K [(1 + 2y) ln(l + 2y) - 2(1 + y) In 2] 



+ 



72 



;i + 4 2/ )(ln(l + 2 2/ )) 2 -4(l+t/)(ln2) 



0(e d 



which is fully analytic in y: 

oo 

F(y) = yJ2 c kV k - 

k=0 

Of course y~ s F(y) and F(x) match in the e-expansion up to second order. 



(32) 



(33) 



4. Parametric form of the equation of state 

The full range of the variables r, M (including the active phase for r < 0) can be 
investigated by introducing a parametric representation for the equation of state. For 
equilibrium systems such a representation was first suggested by Schofield et al. [^6) and 
Josephson |27j and later derived in an e-expansion by Brezin et al. |2l| . To apply this 
method to directed percolation we write r and M as 

br = R(l-6) M = R^O/Oo. (34) 

Inserting ( |34"D with do = 2 into (^5|) one finds 

aq = {R f3 /2) 5 H(6) (35) 
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with the simple scaling function 

H(9) = 9(2-9) + 0(e 3 ). (36) 

The parameter range _R>0,0<6'<2is required to describe the whole phase 
diagram around the critical point. The case r < 0, q — > +0 corresponds to the limit 
# — > 2 (from below). (See figure |l|.) 



Figure 1. Sketch of the parametric representation (p4|). The straight lines correspond 
to a variation of 8 at fixed R. 



M 




T 



As a simple application of the parametric representation we briefly discuss the 
susceptibility 

dM 



X 



dq 



(37) 



which satisfies a power law for q — > 0, 

X = X±\r\~ 1 (38) 

where x+ an d X- correspond to the cases r > and r < 0, respectively and 7 = (3(5—1). 
To second order in e the susceptibility can be written in terms of the parameters 
as 



X 



a2 s-i R -i. 



1 



[1-0)9 



(359(2 - 9) + 2(1- 0)2' ( 39 ) 
Therefore the universal amplitude ratio X-/X+ can be expressed to this order by the 
exponent (3 as 



x+ 



2p-l + 0(e 3 



1 



11 53 , 4\ 
288 - 144 ln 3 )e + ° (e 



1 - -(1 + e 0.0676885...). 



(40) 
(41) 
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The parametric equation of state allows us also to express the expansion coefficient 
c in fl33|) to second order in e by 5. For this purpose we write 

>-<Si 



« = 2W S -W(l- 9 ) V=j$p z = 6->H(6) -z=^§; (42) 
which gives to 0(e 2 ) 



F(y) z 



yl-Sl 



y y (l-O)P(S-i)- 
For y ~ 9 — > this becomes 



2 



(43) 



c = 2 2 ~ d + 0(0 = 1 - - In 2 - — (85 + 164 In 2 - 106 In 3) In 2 + 0(0 (44) 

3 864 

= 1 - e 0.231049 - e 2 0.0659639 + 0(e 3 ). 



5. Conclusion 



We have derived a universal equation that describes the order parameter of directed 
percolation systems as a function of the scaling fields near the transition point to second 
order in e = 4 — d. Using a parametric representation for the thermodynamic variables 
our result could be written in a very simple form. This analysis was motivated by the 
success of a parametric equation of state [analogous to fl3"4])-(|3"6D] for equilibrium critical 
phenomena which was shown to be in good agreement with experiments ||26fl . Our hope 
is that our results will help to analyse data obtained by Monte Carlo simulations. It 
would be especially interesting to see how important the 0(e 3 ) corrections to Eq. (j36|) 
are in practice. 
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Appendix 

In this appendix we give some details of the diagrammatic analysis required to compute 
the equation of state. The diagrams contributing to two loop order are given in figure [AT] . 
The gaussian propagator follows from the action Jq in equation (]8|) as 

q 

t 5 — <— C t' = 6(* - exp (-A(f + q 2 )(t-t')) (1) 



10 



where 0(i) denotes the step function and q is the momentum carried by the line. The 
correlator reads 



t > 

and the vertices 



-c t' 



gM 



2(f + q 2 ) 



cxp (-A(f + q 2 )\t-t' 



(2) 



= \g 

represent the interaction (|9|). 

After integration over the internal time arguments the Feynman diagrams are 

Xg 2 M 



(3) 



1 



(b) 



4 Jpf + p 2 e(2 - e 
Xg 4 M r 1 



■\g 2 Mf x ^l 2 



12 h, P (f + k 2 )(f + p 2 )(f+(k + p) 2 ) 12 
\g 4 M f 1 \g A M _ x 



(d) 



8 h, P (f + k 2 ) 2 (3f + k 2 +p 2 + (k + p) 2 ) 
Xg 5 M 2 r 1 



16 7k, P (f + P) 2 (f + p 2 )(f + (k + p) 2 ) 
A# 5 M 2 1 9 r 1 



(4) 

f L - e h (5) 
r 1 ^ (6) 
(7) 



16 3df 7k, P (f + k 2 )(f + p 2 )(f + (k + p) 

where we have used the notation J* k • • • = (27r)~ d / d d fc • • •. 
The e-expansions of the integrals I\ and I2 read 



h 



Xg 5 M 2 l-e fT 

— f^h 

16 3 



k, P (l + P)(l+p 2 )(l + (k + p) 2 ) 
6G 2 



1 + | + j(7 + /) + 0(e 3 ) 



with 



and 



1 ln(x(l — x)) 

■ da > n V 

1 — x(l — x) 



-2.3439072.... 



(8) 



(9) 



k, P (l + P) 2 (3 + A; 2 + p 2 + (k + p) 1 



3G 2 
2e 



l + -(3-ln3) + 0(e 2 ) 



(10) 



The finite parts of the integrals are required to calculate the coefficient K in the scaling 
function (26) and the normalization constants a and b (28, 29). 
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